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We present a search for charged massive long-lived particles (CMLLPs) that are pair produced in 
pp coUisions at y/s = 1.96 TeV collected by the DO experiment at the Fermilab Tevatron collider. 
Our result is a combination of two searches where either one or both CMLLPs are reconstructed 
in the detector. We select events with muon-like particles that have both speed and ionization 
energy loss (AE/dx) different from muons produced in pp collisions. In the absence of evidence for 
CMLLPs corresponding to 6.3 fb~^ of integrated luminosity, we set limits on the CMLLP masses in 
several supersymmetric (SUSY) models, excluding masses below 278 GeV for long-lived gaugino-like 
charginos, and masses below 244 GeV for long-lived higgsino-like charginos at the 95% C.L. We also 
set limits on the cross section for pair production of long-lived scalar tau leptons that range from 
0.04 pb to 0.008 pb for scalar tau lepton masses of 100 to 300 GeV. 



I. INTRODUCTION 

Several extensions of the standard model (SM) includ- 
ing some SUSY models predict the existence of massive 
long- hved particles (MLLP) [H. Their existence could ex- 
plain the origin of dark matter. Primordial lithium abun- 
dance is not described by the current model of big bang 
nucleosynthesis, but it can be satisfactorily explained by 
the existence of a MLLP that decays during or after big 
bang nucleosynthesis |^ . MLLPs could have color or elec- 
tric charge. They appear as R-hadrons (bound states of 
squarks or gluinos with SM quarks), as sleptons, or as 
charginos. MLLPs are relatively slow moving at the col- 
lision energy of ^/s = 1.96 TeV and for MLLP masses of 
100 GeV or greater considered in this article. Charged 
MLLPs also have large ionization energy loss [AE/Ax) 
due to their slow speeds. These characteristics are differ- 
ent from other particles studied at high energy colliders, 
and thus the identification of such particles is simplified 
by the corresponding small amount of background. We 
therefore search for charged massive long-lived particles 
(CMLLPs) at the Tevatron. 

Searches for CMLLPs were performed previously by 
the DO 041, CDF [1, 0, LEP [|, CMS 0, and AT- 



*with visitors from "Augustana College, Sioux Falls, SD, USA, 
''The University of Liverpool, Liverpool, UK, =UPIITA-IPN, Mex- 
ico City, Mexico, ''DESY, Hamburg, Germany, '^SLAC, Menlo 
Park, CA, USA, ^University College London, London, UK, ^Centro 
de Investigacion en Computacion - IPN, Mexico City, Mexico, 
''ECFM, Universidad Autonoma de Sinaloa, Culiacan, Mexico and 
*Universidade Estadual Paulista, Sao Paulo, Brazil. 



LAS [T^ collaborations. We present limits on masses of 
CMLLPs by combining data from a search for pair pro- 
duced CMLLPs performed with 1.1 fb~^ integrated lumi- 
nosity [5| with an analysis based on 5.2 fb~^ integrated 
luminosity [Tlj . The second analysis includes searches for 
either a pair of CMLLPs or a single CMLLP signature 
in an event. This article provides greater detail on the 
analysis and results published in :11J. 

In this study "long-lived" refers to particles that tra- 
verse the entire detector before decaying. Although cos- 
mological observations place severe limits on stable mas- 
sive particles [J [l2| , these limits do not rule out the par- 
ticles predicted by models studied here. We are sensitive 
to CMLLPs with lifetimes longer than 25 ns, with best 
sensitivity for lifetimes longer than 1 /is. 

We compare the results with predictions of sev- 
eral SUSY models. Models with gauge-mediated 
SUSY-breaking (GMSB) always contain a light grav- 
itino/goldstino as the hghtest SUSY particle (LSP) [ij, 
[13. The next-to-lightest SUSY particle (NLSP) could 
be the lightest scalar tau lepton (stau) or the lightest 
neutralino, depending on the model parameters (isl [l6|. 
The GMSB parameters assumed in this paper make the 
stau lepton the NLSP. If stau lepton decays to the grav- 
itino/goldstino are suppressed (the effective coupling to 
the gravitino/goldstino is a free parameter in the model), 
then the stau lepton can live long enough to escape the 
detector and be a candidate CMLLP [l3, [13 • 

Long-lived charginos can occur in models with anomaly 
mediated SUSY breaking and in SUSY models that do 
not have gaugino mass unification, provided the differ- 
ence between the masses of the lightest chargino and 
the lightest neutralino is less than approximately 150 



4 



MeV [T^ [l^l- The chargino can be mostly higgsino or 
mostly gaugino. We treat these two cases separately. 
The analysis strategy is the same as that for the stau 
lepton search. 

In addition, some SUSY models predict long-lived top 
squark NLSPs that hadronize into mesons and baryons 
with long enough lifetimes to be CMLLP candidates [2l[ . 
Hidden valley models predict scenarios where the top 
squark acts like the LSP and has a long lifetime [12, l23j |. 
In these models the top squark forms hadrons that are 
CMLLP candidates. Any SUSY scenario where the top 
squark is the lightest colored SUSY particle can have 
a hadron formed with a top squark that is a CMLLP. 
Colored CMLLPs will hadronize and experience charge 
exchange during nuclear interactions. This effect is taken 
into account in the analyses reported here. 

A brief description of the DO detector is given in Sec. II, 
which is followed in Sec. Ill by a description of the trigger 
and the data used. Section IV describes the theory and 
the signal generation. Section V presents the strategies 
and techniques used in these analyses. Section VI de- 
scribes the search for pairs of CMLLPs and Sec. VII the 
search for single CMLLPs with an integrated luminosity 
of 5.2 fb"^ "Pair" and "single" refer to the number of 
detected particles. In the models we consider, CMLLPs 
are always produced in pairs. Section VIII summarizes 
the earlier search with an integrated luminosity of 1.1 
fb^^. The combined results are presented in Sec. IX. 
Section X summarizes this study. 



II. DETECTOR 

Figure[T]shows the details of the DO detector [2^ which 
consists of three primary systems: a central tracking sys- 
tem, calorimeters, and a muon spectrometer. The polar 
angle 9 is defined such that 9 = is the +z direction, 
which is the direction of the proton beam. The azimuthal 
angle </> is defined such that (p — lies along the horizontal 
+x axis, pointing outwards from the center of the Teva- 
tron ring and (f) = tt/2 in the +y direction. The pseudo- 
rapidity of a particle is defined as ry = — In [tan (9/2)]. 

The silicon microstrip tracker (SMT) is the innermost 
part of the tracking system and has a six-barrel longitu- 
dinal structure, where each barrel consists of a set of four 
layers arranged axially around the beampipe to measure 
the r-(j) coordinates of charged particles. A new layer of 
SMT sensors was installed near the beampipe in 2006. 
The data recorded before this addition are referred to as 
Run Ila and the subsequent data are referred to as Run 
lib. Twelve radial disks, interspersed between the barrel 
segments, provide position measurement in the r-z and 
r-cj) planes. The SMT provides a spatial resolution of ap- 
proximately 10 /im in r-(j) and approximately 100 /im in 
r-z and covers a pseudorapidity range \r]\ < 3. The SMT 
is also used to measure ionization energy loss (dE/dx) 
of tracks. The central fiber tracker (CFT) surrounds the 
SMT and consists of eight concentric carbon fiber bar- 



rels holding doublet layers of scintillating fibers (one ax- 
ial and one small-angle stereo layer) with the outermost 
barrel covering the region I77I < 1.7. A superconduct- 
ing solcnoidal magnet surrounds the CFT and provides 
a uniform 1.9 T axial magnetic field. 

A liquid argon/uranium calorimeter measures both 
electromagnetic and hadronic energy and is housed in 
three cryostats, with the central calorimeter covering the 
region \ri\ < 1.1 and two end calorimeters covering the re- 
gion 1.5 < 1 77 1 < 4.2. The calorimeter is made of pseudo 
projective towers consisting of an absorber plate and a 
signal board. Liquid-argon, the active material of the 
calorimeter, fills the gap. There are about 10 hadronic 
interaction lengths in the calorimeter at rj = 0. 

The muon system is the outermost par t of the DO de- 
tector and covers the region \r]\ < 2 [2^. It comprises 
drift tubes and scintillation counters arranged in three 
layers (A, B, and C). Between layers A and B, there is 
magnetized steel (6 interaction lengths at ry = 0) gener- 
ating a 1.8 T toroidal field. In the central layers (|7/| < 1) 
multiwire proportional drift tubes (PDT) and in the for- 
ward layers (1 < \ri\ < 2), mini drift tubes (MDT) are 
used for tracking. Scintillation counters covering the re- 
gion (1 77! < 2) are used for triggering on muons. 

The PDTs are typically 2.8 x 5.6 m^, with cells that 
are 10 cm in diameter. Typical chambers are built of 
three or four layers of 24 cell wide planes. Each cell has 
an anode wire at its center. Vernier cathode pads are 
located on both sides of the wires to provide information 
on the hit position along the wire. The chambers are 
filled with a gas mixture of 84% argon, 8% CF4 , and 8% 
CH4 with a drift velocity of approximately 10 cm//Lts. 

Scintillation counters are installed on the top, the sides 
and the bottom of the outer layers of the central muon 
PDTs. They provide a fast signal to associate a muon in 
a PDT with the appropriate bunch crossing and hence 
are used in muon triggers. They also help to discrim- 
inate against the cosmic ray background and to reject 
out-of-time particles scattered from accelerator and de- 
tector components at high rj. The time resolution is ap- 
proximately 2 ns for A-layer counters and approximately 
4 ns for B, and C-layer counters. Detection efficiency is 
close to 100% in all counters. 

In the forward region, MDTs with a drift time of ^ 90 
ns provide good coordinate resolution of less than 1 mm, 
radiation hardness, high segmentation, and low occu- 
pancy. Each MDT layer is divided into octants. An MDT 
consists of eight cells, each with a 9.4 x 9.4 mm^ internal 
cross section and uses a fast gas mixture of CF4/CII4 
(90%:10%). There are 4214 scintillation counters in the 
forward region, arranged in three layers (A, B, and C). 
The segmentation is 4.5° x 0.12 (0.07) in ^ x 77 for the 
first nine inner (last three) rows of counters. The scin- 
tillation counters are 1.3 cm thick with various cross sec- 
tions ranging from 60 x 106 cm^ to 17 x 24 cm^. The 
time resolution is approximately 2 ns and the detection 
efficiency is above 99.9%. The CMLLPs considered in 
this analysis would be identified as muons in the DO de- 
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FIG. 1: (color online) Diagram of the DO detector showing the locations (blue crosses) where a top squark hadron must be 
measured as charged to be selected as a CMLLP candidate. 



tector as they penetrate the material of the calorimeter 
and the toroid and leave hits in the muon system. An 
accurate measurement of the time of flight (TOF) of a 
charged particle reaching a scintillation counter is ob- 
tained from the position of the counter and the recorded 
time of the hit. Particle tracks are reconstructed in the 
muon system using hits from scintillation counters and 
drift tubes. A muon candidate is qualified as a good 
muon if it has hits in scintillation counter layers A and 
either B or C, and multiple drift tube hits in different 
detector layers. These local muons reconstructed by the 
muon spectrometer are then matched to charged particle 
tracks in the central tracking system originating at the 
pp interaction vertex. The muon candidate is rejected if 
no match is found. Otherwise, the measurement of the 
momentum component transverse to the beam line (pt) 
of the muons is taken from the parameters of the cen- 
tral track. To discriminate between muons produced in 
hadronic decays (which tend to be surrounded by other 
charged particles and calorimeter energy deposits) and 
isolated muons, two different isolation quantities are cal- 
culated. Track isolation is the sum of the px of all other 
tracks in the central tracking system in a cone of radius 
AB. = ^(Ary)2 + (A(/))2 < 0.5 around the central track 
matched to the muon. Calorimeter isolation is the sum 
of all energy deposits in the calorimeter in an annulus of 
0.1 < AR < 0.4 around the muon trajectory. Through- 
out this article, isolated muons will be those with track 
and calorimeter isolation less than 2.5 GeV each. 



III. TRIGGER 

The DO trigger system is designed with three distinct 
levels with each succeeding level examining fewer events 



in greater detail so that the final trigger rate is low 
enough for the data to be recorded without causing ex- 
cessive dead time. The first stage (LI) comprising a col- 
lection of hardware elements that selects events based on 
features such as momentum, energy and particle type, 
provides an accept rate of about 2 kHz. In the second 
stage (L2) microprocessors associated with specific sub- 
detectors provide information to a global processor to 
construct a trigger decision based on individual objects 
as well as object correlations. The L2 system has an ac- 
cept rate of approximately 1 kHz. Candidates accepted 
by LI and L2 are sent to the third level (L3) of the trig- 
ger system where the data is processed by algorithms on 
a computing farm to reduce the rate to about 200 Hz. 
Events that pass all three trigger levels are recorded for 
offline reconstruction. 

As described earlier, the CMLLPs in this study have 
muon-like signatures with regards to their penetration 
characteristics in the detector. The px threshold of the 
muon triggering the event varies from 8 GeV to 13 GeV. 
The earlier analysis performed with an integrated lumi- 
nosity of 1.1 fb^^ from Run Ha searched for a pair of 
muons in an event using a triggering condition that re- 
quires at least two muons to be present in the event 
(dimuon trigger). In the recent analysis with an inte- 
grated luminosity of 5.2 fb~^ of Run lib data, the search 
is expanded to include events with only one CMLLP can- 
didate. Hence, triggering conditions requiring the pres- 
ence of at least one muon (single muon triggers) are used. 

Since these triggers are designed for muons traveling 
close to the speed of light (/? ~1), there is a loss of ac- 
ceptance for CMLLP candidates with /? < 1. The muon 
triggers impose a time window (trigger gate) on the muon 
scintillation counter hits. The trigger gate opens 15-30 
ns before particles from a collision traveling at the speed 
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of light reach a particular muon layer and closes 15-40 ns 
after that time. The effect of the trigger gate on events 
with a pair of CMLLPs is more pronounced in di-muon 
triggers, where both slow moving particles must arrive 
within the trigger gate; see Fig. [2ja). The muon trig- 
gers are simulated by applying trigger efficiencies mea- 
sured in Z — > fifi decays in data using a tag and probe 
method. Trigger efficiencies for CMLLPs are calculated 
using Monte Carlo (MC) simulations for different CM- 
LLP masses. The overall selection efficiency, which is a 
product of the trigger gate efficiency and the efficiency for 
single muon triggers, for slow moving CMLLPs is higher 
for single muon triggers than that for dimuon triggers. 
We therefore use single muon triggers in the searches for 
a pair of CMLLPs as well as for a single CMLLP. The 
overall selection efficiency, including the efficiency of the 
trigger gate, for a single CMLLP is shown in Fig. [2][b). 

A second time window (readout gate) is imposed on the 
muon signals during digitization. This gate opens 15-30 
ns before particles from a collision traveling at the speed 
of light reach a particular muon layer and closes 70-90 ns 
after. In the search for a pair of CMLLPs both particles 
must be within the readout gate for the information to be 
recorded even though only one muon needs to be within 
the trigger gate for the trigger to be satisfied. 



IV. MODELS AND SIGNAL GENERATION 

Signal samples with direct production of a pair of CM- 
LLPs are simulated using pythia 6.409 [26]. Data events 
collected from random beam crossings are overlaid on 
simulated events to simulate additional interactions and 
detector noise. Production of CMLLPs through cascade 
decays from heavier new particles (such as squarks) is 
model dependent and has not been considered here. A 
GMSB model with a long-lived stau lepton NLSP, model 
line D in Ref. [27j, is used to generate stau lepton pairs. 
The model parameters are given in Table ID The minimal 
supersymmetric standard model is used for generating 
long-lived gaugino-like charginos, higgsino-like charginos, 
and top squarks. The corresponding model parameters 
are given in Table |lll where Mi , M2 , and M3 are the 
mass parameters for U(l), SU(2), and SU(3) gauginos 
respectively, tan/3 is the ratio of the vacuum expecta- 
tion values of the two Higgs doublets, and fi is the cor- 
responding mass parameter. Long-lived top squarks are 
generated with pythia and hadronized by linking with 
an algorithm external to pythia. This algorithm [28*1 
is applicable to any SUSY model that features a long- 
lived top squark. A set of 50,000 events is generated 
for each model and each mass point. A GEANT-based 
detector simulation models the detector response for the 
MC samples [2^. We have modified GEANT to treat our 
long-lived signal particles as heavy muons for purposes of 
tracking and estimating the dE/dx of signal particles in 
the detector. Therefore, CMLLPs in MC samples have 
muon-like lifetimes. The detector geometry is described 



in detail in the simulation, which uses information on 
the position of the scintillation counters to evaluate the 
timing information of the hits in the counters. This in- 
formation is used to calculate the TOF of the muons. 
Simulation of muon timing in the standard simulation 
software is corrected using information from data as de- 
scribed in Sec. IV Al After the simulation of the detector 
response, a simulation of the electronics and digitization 
is performed. The simulated samples are then passed 
through the same reconstruction software that is used to 
reconstruct data. 

Theoretical values of masses and couplings for dif- 
ferent types of CMLLPs are calculated using Soft- 
SUSY [sol- This information is provided as input to 
PROSPINO 'Sl'l for the calculation of production cross 
sections and their uncertainties for different types and 
masses of CMLLPs. 



V. ANALYSIS STRATEGY AND TECHNIQUES, 
AND SELECTION VARIABLES 

With respect to their production at the primary vertex 
and penetration characteristics, CMLLPs are similar to 
prompt muons produced in pp collisions, but they travel 
at /? « 0.6 to 0.8. The momenta of CMLLPs are dis- 
tinctly higher than that of prompt muons despite their 
lower /3, as shown in Fig. [3] In this figure we compare the 
highest and the second-highest pr CMLLPs (simulated 
gaugino-like charginos of masses 100 GeV, 200 GeV, and 
300 GeV) in an event with the highest and the second- 
highest pt muons from Z ^ decays in data and MC 
events. Events are required to have two isolated muons 
with Pt > 20 GeV. The muons are required to be in- 
consistent with cosmic ray muons. All distributions are 
normalized to the same number of events. For Z — /i/i 
events, we require 70 < M^^ < 110 GeV. The data and 
MC are compared in this mass range to determine the 
corrections to be applied to the MC. The CMLLPs con- 
sidered here have a dE/dx approximately proportional to 
l//3^ which is much higher than the dE/dx of an iden- 
tified muon that is essentially a minimum ionizing par- 
ticle [111. Thus, TOF, which is used to calculate /3, and 
dE/dx of a particle can discriminate between CMLLP 
candidates and muons with /3 ~ 1. 

The Run Ila analysis, which searched for a pair of CM- 
LLPs, used only the TOF of CMLLP candidates to dis- 
tinguish them from prompt muons. The introduction of 
dE/dx measurement into the Run lib analyses allowed 
us to extend the search to events where only one CMLLP 
candidate could be detected. The Run lla dataset was 
used to search for stau leptons, gaugino-like charginos, 
and higgsino-like charginos. The Run lib analyses also 
searched for top squarks. Additional criteria needed for 
the selection of top squark candidates are described in 
detail in Sec. V C. 

Results are presented in this article for stau leptons, 
gaugino-like charginos, and higgsino-like charginos from 
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FIG. 2: (color online) Efficiency for slow-moving gaugino-like charginos of various masses to arrive within the LI muon trigger 
gates. All events in this MC sample contain two gaugino-like charginos. (a) The black and red lines show the efficiencies for a 
pair of charginos, or a single chargino respectively, to be within the trigger gate, (b) Overall efficiency, which is a product of 
the trigger gate efficiency and the efficiency for single muon triggers, for the selection of single charginos. 



Parameter Description Value 

A™ Scale of SUSY breaking 19 to 100 GeV 

Mm Messenger mass scale 2 A™ 

Ni, Number of messenger fields 3 

tan P Ratio of Higgs field vacuum expectation values 15 
sign of /i Sign of Higgsino mass parameter -1-1 



TABLE I: GMSB model parameters for stau lepton production. 



Model /i (GeV) A/i (GeV) Ma (GeV) M-j. {GeV) tan^ Squark Mass (GeV) 



Top squark 


10,000 


100 


200 


500 


15 


800 


Gaugino-like chargino 


10,000 


3M2 


100 to 300 


500 


15 


800 


Higgsino-like chargino 


100 to 300 


100,000 


100,000 


500 


15 


800 



TABLE II: Model parameters for top squark and chargino production. 

10'= 
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FIG. 3; (color online) Distributions of pr of the (a) highest pr and the (b) second-highest pr muon in an event for Z — > fj,^, 
data and MC, and for simulated gaugino-like charginos with masses of 100 GeV, 200 GeV, and 300 GeV. 
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individual analyses as well as from a combination of Run 
Ila and Run lib data with a total of 6.3 fb~^ integrated 
luminosity. To avoid double counting of events, the sam- 
ples used for the two Run lib analyses are constructed 
to be statistically independent. All events that pass the 
selection requirements used for the search for a pair of 
CMLLPs are removed from the data used for the single 
CMLLP search, resulting in an approximately 40% loss 
of signal acceptance for the single CMLLP search. The 
data and the background sample that are used in the 
search for single CMLLPs contain muons that originate 
mostly from the decays of W bosons. The number of 
such events changes by only about 2% due to this veto. 
Furthermore, we show results for searches for single top 
squarks and top squark pair production using 5.2 fb~^ 
integrated luminosity. A combination of the two analy- 
ses does not improve the result as explained in Sec. V C, 
and is not performed here. 



A. Time-of-Flight Measurement 

The TOF of a charged particle reaching the muon sys- 
tem can be calculated from the position of the scintilla- 
tion counter that is hit and the corresponding time as 
has been described in Sec. II. A time offset is determined 
for each scintillation counter along with its associated 
cables and time digitizing electronics using a sample of 
muons from experimental data so that a time reading of 
zero is obtained for particles that originate at the center 
of the detector at the time of a beam-beam interaction 
and travel at the speed of light to the specific scintilla- 
tion counter. These time offsets are imperfect and have 
fluctuated with, for example, the seasonal variation in 
the synchronization with the Tevatron accelerator clock 
±1 ns) as shown in Fig. HJ^a). The offsets are cor- 
rected by subtracting the relevant amount of deviation 
from each hit time in each scintillator. FigureH^b) shows 
the mean of the time distribution of hits versus time af- 
ter the offset correction. The time period for averaging 
is the duration for which p and p beams are circulated in 
the Tevatron accelerator after injection. This period is 
typically 12-24 hours. 

It is observed that the MC simulation gives narrower 
time distributions than what is observed in data. There- 
fore, the TOF associated with a muon hit in MC is 
smeared to reflect the resolution of the time distribution 
of muons from Z boson decays in data. The amount of 
smearing depends on the location of the muon detector 
because of the differing sizes of the scintillation counters. 
Figures [5][7] show the time distributions in selected re- 
gions for layers A, B, and C in data and MC. There is 
some mismodeling at early times in the central C layer 
(Fig.mja)) and central bottom B layer (Fig. [Tljc)) arising 
from the data-driven smearing that has been applied to 
the hit times of the muons in MC. This mismodeling is 
not in our signal region which is at large times and has 
a negligible effect on the results. 



The value of /3 for a muon-like track is determined from 
a weighted average of the speeds Pi determined using 
times corresponding to individual scintillation counter 
hits on the track, the weights being the inverse of 
the squares of experimentally determined uncertainties. 
Figure [5] shows the distribution of the highest and 
the second-highest px muons in data and in simulated 
gaugino-like chargino events. 

We define the speed significance, which incorporates /3 
with its uncertainty icr{^) ) as. 



speed significance — 



(1) 



Figure [9] shows speed significance distributions of the 
highest and the second-highest pr muons in the event. 
Prompt muons will have speed significance near zero 
whereas CMLLPs will have positive speed significance. 
In the Run Ila analysis we use the product of the speed 
significances of the two muons, which will be positive for 
a pair of CMLLPs, as an additional criteria to separate 
signal events from background events. We also calculate 
the variable. 



speed x^/dof 



TV- 



m - A 



(2) 



where N is the number of scintillation counter hits as- 
sociated with a muon track, /?i is the speed of the track 
corresponding to the hit in the i*'* scintillation counter 
and (/?) is the weighted average of the speeds calculated 
for all scintillation counter hits on the track. For some 
tracks the measurement /3i from the hit on a particu- 
lar layer is inconsistent with the /3j measurements from 
other layers causing a large contribution to the value of 
the speed x^/dof. To identify and remove this hit from 
the speed x^/dof calculation, hits on a muon track are 
removed one at a time and the speed x^/dof is recalcu- 
lated with the rest of the hits. The set of hits with the 
lowest value of speed x^/dof is used to recalculate the 
speed provided it satisfies the qualities of a good muon 
candidate (described in Sec. II). Distributions of speed 
X^/dof for the highest px muon in data and in simulated 
gaugino-like chargino events are shown in Fig. [TOja) . 

The speed asymmetry, which is useful in the search for 
a CMLLP pair, can be defined for the two highest pT 
muons in the event as: 



speed asymmetry 



/3i+/32 



(3) 



The speed asymmetry is near zero for both signal events 
and events containing two well measured muons, but will 
be large for events where one of the particles has a speed 
that has been mismeasured. Speed asymmetry distri- 
butions for data and for simulated gaugino-like chargino 
events are shown in Fig. [TU] (b) . 

We have observed that the disagreements in (P) distri- 
butions from data and MC at small values of (/?) , visible 
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FIG. 4: (color online) Mean hit time as a function of the Tevatron store number in the forward region A-layer for muons from 
the decay of Z bosons in data (a) before correction and (b) after correction. A store is the time period for which p and p beams 
are circulated in the accelerator after injection. This period is typically 12-24 hours for the Tevatron. When the number of 
recorded times from Z — >■ up decays for a given store is large, a Gaussian function is fit to those times and the Gaussian mean 
is used. If the number of recorded times for a given store is small, the median of these times is used. These distributions cover 
the data taking period between June 2006 and March 2010. 
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FIG. 5: (color online) Time distribution for scintillation counter layers A, B, and C in the forward muon system for times from 
Z — !> pp decays for data and MC events. 
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FIG. 6; (color online) Time distributions for scintillation counters in layer C, (a) top and sides, (b) bottom, in the central 
muon system for times from Z — >■ \ni decays for data and MC events. 
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FIG. 7: (color online) Time distributions for scintillation counters in (a) layer A (top, bottom and sides), (b) layer B sides, and 
(c) layer B bottom, in the central muon system for times from Z — )■ \ni decays for data and MC events. 
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Highest Muon's < p > Second Highest IVIuon's < p > 



FIG. 8: (color online) Distributions of (/3) of the (a) highest pr and the (b) second-highest pr muons. The distributions are 
normalized to the same number of events. Background is taken from Z — >■ fifi MC decays. The selection requirements are 
identical to those used in the pair CMLLP analysis, as described in Sec lVIl except that the requirement (/3) < 1 is not imposed. 




Highest p.^ Muon's Speed Significance Second Highest p.^ Muon's Speed Significance 



FIG. 9: (color online) Speed significance of the (a) highest pr and the (b) second-highest pr muons. The distributions are 
normalized to the same number of events. Background is taken from Z — >■ ^/j. MC. The selection requirements are identical to 
those used in the pair CMLLP analysis, as described in Sec IVIl except that the requirement {(3) < 1 is not imposed. 



in Fig. IHl can be removed by applying stringent require- 
ments on either the speed x^/dof or the speed asymmetry 
distribution. The (/?) distribution extends beyond 1 due 
to uncertainty in the measurement of velocity. There is 
some mismodeling for > 1, but this is not in our sig- 
nal region. This disagreement is due to the background 
coming from mismeasured muons, which is characterized 
by large values of speed x^/dof or speed asymmetry, and 
not due to signal-like events, which are characterized by 
small values of speed x^/dof or speed asymmetry. We 
apply a requirement on the speed asymmetry of the two 
candidate CMLLPs in the search for a pair of CMLLPs 
as described in Sec. VI B. In the search for single CM- 
LLPs we apply a requirement on the speed x^/dof of the 
candidate CMLLP as described in Sec VII B. 

We correct the mismodeling in the speed x^/dof and 
the speed asymmetry distributions using a signal-free re- 
gion (as described in Sec. VII B) in data. We compute 



the ratio of the speed x^/dof (or the speed asymmetry) 
distributions in data and Z — >■ /x/x MC for events with 70 
< Mfj,fj, <110 GeV, where the potential signal contribu- 
tion is negligible. The value of this ratio is applied as a 
weight to all simulated events. 



B. The dE/dx Measurement 

The dE/dx of a particle is measured in the SMT. It is 
a calibrated average over SMT clusters and is corrected 
for the path length of the particle in the barrel or the disk 
sensor to reduce the dependence on the incident angle of 
the particle. The calculation of dE/dx excludes SMT 
clusters with the highest 20% of dE/dx values in order 
to minimize the contribution from Landau tails. 

The average value of dE/dx is observed to decrease 
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FIG. 10: (color online) (a) Distribution of speed x^/dof of the highest pr muon in an event, and (b) speed asymmetry 
distribution of the highest and the second-highest pr muons for Z — > data, MC and signal (gaugino-like chargino with 
masses of 100 GeV, 200 GeV, and 300 GeV). The distributions are normalized to the same number of events. All events are 
required to have two isolated muons with pr >20 GeV. For Z — >■ events, we require 70 < M^^ < 110 GeV. The muons are 
required to not be consistent with cosmic ray muons. 




FIG. 11: (color online) Distributions of the adjusted dE/dx (described in Sec. V B) of the (a) highest pr and the (b) second- 
highest pt muons in data and the simulated chargino events. The distributions are normalized to the same number of events. 
The Background is taken from Z — >■ fifi MC. The selection requirements are identical to those used in the pair CMLLP analysis, 
as described in Sec lVIl 



with increasing integrated luminosity [33[ due to radia- 



tion damage to the silicon sensors. To correct for this 
effect, each dE/dx measurement is divided by the mean 
di?/da;.at dE/dx — 1, to facilitate the comparison be- 
tween data and MC. The recalibrated dE/dx distribution 
is referred to as the "adjusted dE/dx'\ Figure [TT] shows 
that the adjusted dE/dx distribution of muons from Z 
decays is well separated from that of candidate CMLLPs. 
Since the adjusted dE/dx distribution of muons from Z 
decays in data does not quite match the adjusted dE/dx 
of Z — > /Lt/i MC events, a Gaussian smearing is applied to 
the adjusted dE/dx in MC to better describe the data. 
The precision of a particle's adjusted dE/dx depends on 
the number of SMT clusters used in its calculation. A 
new variable, dE/dx significance, is defined to estimate 



this dependence. If Nc is the number of SMT clusters on 
a track, and if the spread in the adjusted dE/dx distri- 
bution (which is a function of Nc) is a{dE/dx)N^, then 
the dE/dx significance is: 



dE/dx significance = 



d£;/da; - 1 
a{dE/dx)N, 



(4) 



Figure [12] shows a distribution of the dE/dx significance 
for the highest px and the second-highest pr muons 
in data and background, and for simulated gaugino-like 
charginos of 100-300 GeV masses. The two data points 
at high values of dE/dx significance have (/3) close to 
0.98 and are removed by the selection criteria described 
in Sec. VI B. 



13 




Highest Muon's dE/dx Significance Second IHighest IVIuon's dE/dx Significance 



FIG. 12: (color online) Distributions of dE/dx significance for the (a) highest pr and the (b) second-highest pr muons. The 
distributions axe normalized to the same number of events. Background is taken from Z — > nfi MC. The selection requirements 
are identical to those used in the pair CMLLP analysis, as described in Sec lVIl 



C. Detection of Top Squarks 

About 60% of top or anti-top squark hadrons will be 
charged after initial hadronization j34| . A top squark 
hadron passing through matter can exchange light quarks 
through nuclear interactions, changing the charge of the 
hadron. Because the detector has more quarks than an- 
tiquarks, after many nuclear interactions most of the top 
squark mesons become baryons, but anti-top squark anti- 
baryons become anti-mesons. The charge of a top squark 
baryon can be 0, -f 1 or +2; the fractions are model de- 
pendent. We have assumed 2/3 of the stop baryons to be 
charged after undergoing many interactions in the detec- 
tor material. Similarly, the charge of an anti-top squark 
anti-mesons can be or —1. We assume anti-top squark 
anti-mesons to have a probability of 1/2 of being charged 
after passing through the detector material |35l - l37| . We 
also provide results assuming the top and anti-top squark 
hadrons do not flip charge at all for reference. In this case 
we only include a factor of 60% for the initial hadroniza- 
tion. 

As can be seen in Fig. [1] the top squark or the anti- 
top squark hadron must be measured as charged at three 
locations while passing through the DO detector to be re- 
constructed: after hadronization, after the calorimeter, 
and after the muon toroid. Both the calorimeter and the 
muon toroid contain enough material (10 and 6 inter- 
action lengths respectively) to ensure that a top squark 
or anti-top squark hadron will undergo a large number 
of interactions when passing through them, randomizing 
its charge. Therefore, the probability of a top squark 
hadron to be charged at all three locations is 0.6 (at pro- 
duction) X 0.67 (at the end of the calorimeter) x 0.67 
(at end of the muon toroid) = 0.27. Likewise, the prob- 
ability of an anti-top squark hadron to be charged at all 
three locations is 0.6 x 0.5 x 0.5 = 0.15. The probabil- 
ity for a pair of top squark and anti-top squark hadrons 



both to be charged in all three locations is 0.27 x 0.15 
= 0.04. The probability of at least one of the two be- 
ing charged in all three locations is 0.27 x (1 — 0.15) 
+ 0.15 X (1 - 0.27) + 0.27 x 0.15 = 0.38. We apply 
high enough pT cuts on the reconstructed tracks so that 
the selected tracks have small curvatures and even the 
tracks with +2 charge are reconstructed and matched. 
For both the pair and the single CMLLP searches, these 
estimates of charge survival probabilities are applied as 
additional factors when the top squark MC is normalized 
to the expected number of events. 



VI. SEARCH FOR EVENTS WITH A PAIR OF 
CMLLPS 

With the selection variables described in Sec. V to pro- 
vide discrimination of CMLLP signal over background, 
we describe below the selection criteria for the search for 
a pair of CMLLPs with 5.2 fb"^ of Run lib integrated 
luminosity. 



A. Background Sample 

To model the background, one million events contain- 
ing muons from decays of Z bosons were simulated us- 
ing PYTHIA. Although a, Z fi^ MC sample is used 
to model the background, we do not assume that the 
background is only from Z — > /x/i decays. Any source 
of muons that are not measured correctly contributes to 
the background. Since the px of the background muons 
have a wide range, we minimize the dependence of the 
analysis on the px of the selected tracks. Other than the 
initial requirement on the pt of the tracks, the selection 
criteria do not depend any further on the px of the candi- 
date CMLLPs. The important variables in this analysis. 
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p and normalized dE /Ax, are largely independent of the 
Pt of the muon. The simulation has been tuned so that 
the /3 and the d£'/da; are well modeled as described in 
Sec. V A and B. 



B. Event Selection 

Events containing a pair of CMLLP candidates are se- 
lected by requiring that they contain at least one muon 
with pT > 20 GeV and exactly two reconstructed, iso- 
lated muons of good quality with |?7| < 2. The other 
selection criteria are: 

• Number of SMT clusters, (Nsmt) > 3. 

• < 1 for the two highest pr muons. 

• asymmetry < 0.35. 

• The highest pT muon has pr > 55 GeV and the 

second-highest pr muon has pr > 50 GeV. 

• Matching (between the track in the muon sys- 
tem and the central tracker) < 100. 

Muons originating from cosmic rays are vetoed by check- 
ing the difference in their arrival times at the scintilla- 
tion counter layers A and C shown in Fig. [rSf aV The 
difference in the TOFs from the center of the detector to 
layer- A of two muons (Fig. ll3f b)) also provides discrimi- 
nation of muons produced in pp collisions from cosmic ray 
muons. The distance of closest approach (DCA) in the 
r-0 plane of each reconstructed muon track to the beam 
line, and the pseudo-acolinearity, Aa = \ A(j) + A9 — 2tt\, 
between the two muon tracks are also used to reject cos- 
mic ray muons. These criteria are as follows: 

• DCA in the r-cj) plane of each muon < 0.2 cm. 

• Difference between A-layer and C-layer times of a 
muon > — 10 ns. 

• Absolute value of the difference in A-layer times 
between the two muons < 10 ns. 

• Aa > 0.05. 

These selection criteria are optimized to produce the best 
expected limits on the masses of CMLLP candidates. 

The background sample is normalized to data events 
that pass the conditions described above and have in- 
variant mass of the two highest px muons in each event 
within 70 < M^^^, < 110 GeV. The contribution of po- 
tential signal in this region is negligible. The number 
of data and background events selected using the above 
conditions are listed in Table IIIIl and the efficiencies for 
signal events for various CMLLP candidates are given in 
Tables irvUVlIl 

A Boosted Decision Tree (BDT), as implemented in 
Ref. [ssj . is used to further discriminate signal events 
from background events. The BDT is trained using the 



expected signal and background distributions from MC, 
modified as described above. Half of the events are used 
to train the BDT while the remaining half are used to 
test the background model and signal response of the 
BDT. The variables used as inputs to the BDT are /3, 
speed significance, normalized d£'/da;, and dE/dx sig- 
nificance for the highest px and the second-highest pt 
muons. The distributions for these variables are shown 
in Figs. [3 H El and [12 where Figs. M and M show the 
distributions without the (/3) < 1 requirement. The cor- 
relation matrices for these variables for a stau signal with 
a mass of 300 GeV and muons from Z boson decays in 
MC are shown in Fig. [T4l a) and (b) , respectively. These 
figures show that there are non-trivial correlations be- 
tween the variables and therefore MVA methods such as 
BDT, which will be able to take the correlations into ac- 
count, are appropriate for this analysis. An example of 
the BDT output distribution is shown in Fig. [T5] for sim- 
ulated gaugino-like charginos of different masses. The 
BDT outputs for simulated stau leptons, top squarks, 
and higgsino-like chargino signals are shown in Figs. 1241 
[2n]in Appendix A. The final selection criteria on the BDT 
output are optimized to yield the best expected cross sec- 
tion limit in the no-signal hypothesis for each mass point, 
for each signal type. 

C. Systematic Uncertainties 

Systematic uncertainties are included in the estimation 
of cross section limits as follows. Each input parame- 
ter to the BDT distribution used to distinguish between 
signal and background is varied within its one standard 
deviation uncertainty and a new BDT distribution is pro- 
duced for each variation for both signal and background 
models. The new BDT distributions are compared to 
the nominal ones and the average of the change in the 
occupancy of the BDT bins satisfying the selection re- 
quirement is taken as the systematic uncertainty due to 
that parameter. This procedure is applied to the sources 
of systematic uncertainties summarized below. These un- 
certainties are used in the limit calculation to model the 
effects of systematics in determining the limits. 

• Uncertainty on the muon momentum scale in simu- 
lated samples is calculated by varying the pr of the 
muon track by ±1 standard deviation as measured 
in data. 

• Time distributions for background and signal are 
modeled using Z — >■ fifi data events. We repeat 
the time smearing with time distributions of muons 
coming from W — )■ fiv decays in data and take the 
difference as the uncertainty. 

• To account for effects of the calibration of individ- 
ual scintillation counters on the LI trigger gate, we 
shift the trigger gate by ±1 ns and calculate the 
resulting change in signal efficiency. 
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FIG. 13: (color online) (a) Difference between the A-layer and C-layer times for a single muon. There are two cosmic ray peaks 
for the two possible directions, away from or towards the pp collision vertex, (b) Absolute value of the difference between the 
A-layer times of the two muons in the event. The times shown in these plots are centered at zero for (/?) = 1 particles. This 
cosmic ray data was collected when there was no p or p beam in the Tevatron collider. Selection requirement on the time 
difference is shown with a blue vertical line. 



Selection Criteria 


Data events 


Background Acceptance 
Z ^ /i/i MC (%) 


Initial muon selection 


231487 


49.6 


N^ = 2 


178204 


15.8 


Trigger matching 


125662 




Trigger probability 




8.85 


Cosmic veto 


106941 


7.57 


LI trigger gate (MC only), Nsmt > 3 


98195 


6.94 


(13) < 1, (13) asymmetry < 0.35 


34376 


2.17 


>55, pt^^ >50 GeV 


709 


0.03 


Matching <100 


702 


0.03 



TABLE III: Selection efficiencies for data and background events before the application of BDT requirements for the search for 
a pair of CMLLPs with Run lib data. Initially, events with at least one isolated muon of good quality and pr > 20 GeV are 
required and a match between a muon and a central track is required based on x^- 



M(stau lepton) in GeV 
Selection criteria 


100 


150 


200 


250 


300 


Initial muon selection 


62.0 


61.4 


62.2 


61.4 


62.0 


N^ = 2 


25.0 


25.2 


25.6 


25.6 


25.5 


Trigger probability 


14.4 


14.6 


14.6 


14.5 


14.4 


Cosmic veto 


12.2 


11.9 


11.7 


11.6 


11.5 


LI trigger gate, Nsmt > 3 


10.5 


10.2 


9.93 


9.50 


8.84 


{/?) <1, iP) asymmetry<0.35 


9.08 


9.44 


9.45 


9.17 


8.62 


pt^} >55, pi;? > 50 GeV 


8.35 


9.30 


9.42 


9.15 


8.60 


Matching < 100 


8.24 


9.15 


9.29 


9.00 


8.45 



TABLE IV: Selection efficiencies (in %) before the application of BDT requirements for a pair of stau leptons in simulated 
events. The initial muon selection requires at least one isolated muon of good quality, matched to a central track with pr > 20 
GeV. 
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M(top squark) in GeV 


100 




150 


200 


250 


300 




350 


400 


Selection criteria 




















Initial muon selection 


50.6 




54.9 


57.9 


59.6 


58.6 




58.3 


57.7 


N^ = 2 


18.1 




20.7 


22.3 


22.9 


21.9 




21.9 


21.2 


Trigger probability 


10.4 




12.0 


12.9 


13.2 


12.5 




12.5 


12.0 


Cosmic veto 


8.67 




9.81 


10.5 


10.5 


9.72 




9.71 


9.31 


LI trigger gate, Nsmt > 3 


7.79 




8.79 


9.26 


8.83 


7.63 




6.76 


5.34 


(P) < 1, (13) asymmetry < 0.35 


7.16 




8.28 


8.93 


8.63 


7.45 




6.65 


5.24 


> 55, > 50 GeV 


5.97 




8.06 


8.85 


8.61 


7.44 




6.64 


5.24 


iviaicning ^ ^iuu 


5.90 




7.96 


8.73 


8.49 


7.34 




6.53 


5.18 


Charge survival efficiency (4%) 


0.24 




0.32 


0.35 


0.34 


0.29 




0.26 


0.21 


TABLE V: Selection efficiencies (in 


%) before the application of BDT requirements for a pair of top squarks in simulated events. 


The initial muon selection requires 


at least one 


isolated muon 


of good quality, matched to a 


central track with pr > 20 GeV. 


The top squark charge survival efficiency is 4% 


















M(gaugino-like chargino) in GeV 




100 




150 


200 




250 




300 


Selection criteria 




















Initial muon selection 




46.8 




49.8 


50.3 




49.1 




49.1 


= 2 




16.2 




17.9 


18.3 




17.6 




17.7 


Trigger probability 




9.32 




10.4 


10.6 




10.3 




10.3 


Cosmic veto 




7.73 




8.46 


8.42 




7.94 




7.77 


LI trigger gate, Nsmt > 3 




6.96 




7.39 


7.11 




6.09 




5.29 


(/3) <1, (13) asymmetry<0.35 




6.36 




6.98 


6.88 




5.95 




5.16 


pt^} >55, plj? >50 GeV 




5.23 




6.67 


6.77 




5.92 




5.14 


Matching <100 




5.16 




6.59 


6.67 




5.85 




5.05 


TABLE VI: Selection efficiencies ( 


in %) before the application of BDT requirements for a 


pair of 


gaugino-like charginos in 


simulated events. The initial muon selection requires 


at least 


one isolated 


muon of good quality, matched to a 


central track 


with pT > 20 GeV. 




















M(Higgsino-Iike chargino) in GeV 




100 




150 


200 




250 




300 


Selection criteria 




















Initial muon selection 




47.5 




50.5 


51.4 




51.1 




49.0 


N^ = 2 




17.1 




18.6 


19.3 




19.0 




18.0 


Trigger probability 




9.78 




10.8 


11.2 




11.1 




11.4 


Cosmic veto 




8.15 




8.87 


8.88 




8.65 




8.04 


LI trigger gate, Nsmt > 3 




7.29 




8.13 


7.56 




6.83 




5.69 


(/3) < 1, (/?) asyrnmctry<0.35 




6.57 




7.70 


7.31 




6.68 




5.59 


p^^ > 55, > 50 GcV 




5.40 




7.39 


7.23 




6.65 




5.56 


Matching < 100 




5.31 




7.31 


7.14 




6.56 




5.49 



TABLE VII: Selection efficiencies (in %) before the application of BDT requirements for a pair of higgsino-like charginos in 
simulated events. The initial muon selection requires at least one isolated muon of good quality, matched to a central track 
with PT > 20 GeV. 
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Correlation Matrix (signal) 
(a) D0 Simulation 



Linear correlation coefficients in % 



Correlation Matrix (background) 
(b) D0 Simulation 
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FIG. 14: (color online) Correlation matrix for different kinematic variables for (a) stau leptons of 300 GeV mass, and (b) 
background for the search of a pair of CMLLPs in the Run lib data. 



• Uncertainty due to the correction to the AEjAx 
measurement of selected tracks to equalize the de- 
grading response due to radiation damage in silicon 
is evaluated by varying the correction factor by its 
±1 standard deviation uncertainty. 

• d-E/dx modeling uncertainty: the Gaussian smear- 
ing function applied to the di?/dx distribution of 
muons in MC is derived separately using Z — >■ /i/i 
and W ^ yaj data events. The difference between 
the BDT distributions obtained with the two meth- 
ods is taken to be the systematic uncertainty due 
to di?/da; modeling. 

• Theoretical values of production cross sections (de- 
scribed in Sec. IV) depend on the choice of parton 
distribution functions (PDF). Their effect is esti- 
mated by using the 40 CTEQ6.1M error PDFs [11] 
for signal and background. The variations from 
each of the error PDF sets and from the renor- 
malization and factorization scale uncertainties are 
added in quadrature. 

The remaining systematic uncertainties, given below, 
are added in quadrature to the uncertainties described 
above to obtain the total systematic uncertainties on sig- 
nal acceptance and background prediction. 

• The uncertainty in muon identification is a com- 
bination of the uncertainties due to selection of a 
muon (1.2%), central track reconstruction (1.4%), 
and isolation of the muon (0.9%). 

• To determine the uncertainty in background nor- 
malization, the mass window for the control re- 



gion is changed from 70 < M^^ < 110 GeV to 60 
< A/^^ < 120 GeV and 80 < M^^ < 100 GeV. 

• The MDTs have an asymmetric timing gate with a 
total length of 94 ns, which is not modeled in the 
MC. The signal from the earliest muon arrives at 
the MDTs within 74 ns of the beam crossing from 
Z ^ data. At most 1.2% of the GMLLPs that 
we consider will be in the forward muon system 
and will arrive at the MDTs more than 20 ns af- 
ter a prompt muon. We have therefore introduced 
an additional 1.2% uncertainty on the signal accep- 
tance. 

• Uncertainty due to the speed asymmetry correction 
(described in Sec. IV A[) is estimated by the change 
in signal acceptance with and without the correc- 
tion. This uncertainty is found to vary between 
1% and 10% depending on the masses and types of 
CMLLPs. The value of this uncertainty is 3.6% for 
the background sample. 

An uncertainty of 6.1% 33] on the integrated luminos- 
ity is applied to the signal efficiency. Systematic uncer- 
tainties for signal and background samples are listed in 
the second columns of Tables IXVlIl and IXVIIII The to- 
tal systematic uncertainty for the background estimate is 
18.2% and that for the signal acceptance is 11.2-15.2%, 
depending on the signal model and mass of the CMLLP. 

D. Results 

The signal acceptance, the background prediction, and 
the observed number of events after the BDT require- 
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FIG. 15: (color online) BDT-output distributions for simulated gaugino-like chargino masses of 100-300 GeV in 50 GeV steps 
for the search for a CMLLP pair with the Run lib data. Distributions are normalized to the expected number of events. 
Selection requirement on the BDT value is shown with a green vertical line. Note that a different BDT was constructed for 
each mass and the BDT selection requirement optimized separately for each mass. 



ment are shown in Table IVIIII These numbers are used 
as inputs to a modified Bayesian method [13] for cal- 
culating the limits on production cross sections at 95% 
C.L. Theoretical values of the production cross sections 
(described in Sec. IV) and observed and expected val- 
ues of limits on the production cross sections of various 
CMLLPs are shown in Table |IX] and Fig. HH The lower 
mass limits that are obtained from the cross section lim- 
its are 189 GeV for top squarks, 250 GeV for gaugino-like 
charginos, and 204 GeV for higgsino-like charginos. The 
limit on the mass of top squarks would increase to 280 
GeV if we would include only the effects of their initial 
hadronization. If the intersection point of the —1 (+1) 
standard deviation band with the NLO cross section is 
used, then the mass limits shift down (up) by ~1 GeV 
for charginos and by ^20 GeV for top squarks. 



VII. SEARCH FOR EVENTS WITH A SINGLE 
CMLLP 

The following section describes the search for a sin- 
gle CMLLP in 5.2 fb~^ of i nteg rated luminosity. More 
details can be found in Ref. [la. 



A. Background Sample 

The dominant background in the single CMLLP search 
is muons from the decays of W bosons, which is modeled 
with data. To define independent data and background 
samples, we select events using the transverse mass, My, 
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TABLE VIII: Signal acceptance, number of predicted background events, and number of observed events in the search for a 
pair of CMLLPs with Run lib data after each BDT selection. The error is a sum in quadrature of statistical and systematic 
errors. 



of the W boson given by 

Mt - ^(PT^ ± ^t)2 - (p. ± fi.r - (Py + f^y?- (5) 

Here is the transverse momentum of the muon and 
f^T is the total unbalanced momentum transverse to the 
beamline as measured in the calorimeter and corrected 
for the muons. Events with Mr < 200 GeV and (/?) < 1 
are selected for the background sample and events with 
Mt > 200 GeV and {(3) < 1 constitute the search sam- 
ple. Figure [17] shows a distribution of Mt for single 
muon events from data and for higgsino-like chargino MC 
events. 



B. Event Selection 

The criteria to select events with one or more CM- 
LLPs are similar to those used in the search for a pair of 
CMLLPs (Sec. VI B). Events satisfying a suite of single 
muon triggers are required to contain an isolated muon 



of good quality (described in Sec. VI B) within |?7| < 1.6. 
The muon must originate at the pp interaction vertex and 
must satisfy the following criteria: 

• PT > 60 GeV. 

• speed x^/dof < 2. 

If there is more than one such muon in the event, only 
the highest pt muon is considered as the CMLLP can- 
didate. To ensure that selected muons do not originate 
from cosmic rays, we require the DC A in the r-cj) plane of 
the selected muon track to the beam line to be less than 
0.2 cm and the difference between A-layer and C-layer 
times to be > — 10 ns. 

If there is a second muon passing all the selection cri- 
teria, conditions to remove the cosmic ray events are the 
same as in the search for CMLLP pairs (described in 
Sec. VI B). The selection criteria and the corresponding 
efficiencies for number of events in data and the CMLLP 
signals are given in Tables |X] - IXIVI 

Since the background is modeled using data, it is nec- 
essary to normalize it to data in a signal-free region. 
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TABLE IX: NLO cross section and 95% C.L. limits (ergs) for the search for a pair of CMLLPs in the Run lib data. The top 
squark charge survival efficiency is 4% (see text). 



Selection Criteria 



Number of Events 



Initial Muon Selection 
Isolated muon 
Trigger matching 

Cosmic veto 

NsMT > 3 

PT > 60 GeV 
{13} < 1, Speed xVdof < 2 
Matching < 100 



345173 
141008 
79983 
79096 
57532 
56466 
27876 
27742 



Background (Mt < 200 GeV) 
22368 



Data (Mt > 200 GeV) 
5374 



TABLE X: Selection efficiencies for data in the search for single CMLLPs before applying the BDT selection requirements. 
Initially, events with at least one isolated muon of good quality and pr > 20 GeV are required and a match between a muon 
and a central track is required based on x^. 
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FIG. 16: (color online) 95% C.L. limits on production cross sections of a pair of stau leptons, top squarks, gaugino-like charginos, 
and higgsino-like charginos as a function of their masses from the search for a pair of CMLLPs with Run lib data. "CF" is the 
scenario without charge flipping. ±1 SD and ±2 SD are the 1 and 2 standard deviation bands respectively around the expected 
limit curves. 



M(stau lepton) in GeV 
Selection criteria 


100 


150 


200 


250 


300 


Initial muon selection 


68.9 


67.9 


68.9 


68.1 


68.5 


Isolated muon 


59.1 


59.5 


60.6 


60.2 


60.8 


Trigger probability 


31.5 


31.2 


31.4 


31.0 


31.0 


Cosmic veto 


27.9 


26.2 


25.3 


27.3 


26.1 


LI trigger gate, Nsmt > 3 


22.2 


19.6 


18.0 


18.0 


15.6 


PT > 60 GeV 


20.3 


19.4 


18.0 


18.0 


15.6 


< 1, Speed xVdof < 2 


16.9 


16.5 


15.6 


15.7 


13.7 


Matching < 100 


16.8 


16.4 


15.5 


15.6 


13.6 


Mt > 200 GeV 


14.4 


15.6 


15.2 


15.4 


13.6 



TABLE XI: Selection efficiencies (in %) for a single stau lepton in MC events. Initially, events with at least one isolated muon 
of good quality and pr > 20 GeV are required and a match between a muon and a central track is required based on x^- 
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• 

M(top squark) in GeV 
Selection criteria 


1 nn 
iUU 


iOU 


zUU 


ZOU 


Qnn 
oUU 


Q pen 


4UU 


Initial muon selection 


56.9 


62.3 


65.6 


56.6 


67.2 


66.9 


66.5 


Isolated muon 


46.0 


52.0 


55.8 


47.0 


57.3 


57.3 


56.6 


Trigger probability 


24.9 


27.6 


29.3 


25.3 


29.4 


29.3 


28.8 


Cosmic veto 


21.1 


23.3 


26.2 


24.3 


24.6 


23.3 


24.6 


LI trigger gate, Nsmt > 3 


15.7 


16.8 


18.0 


15.3 


14.3 


11.8 


10.6 


PT > 60 GeV 


13.0 


16.3 


18.0 


15.3 


14.2 


11.8 


10.6 


{/?) < 1, Speed xVdof < 2 


11.1 


14.3 


15.7 


13.5 


12.5 


10.4 


9.3 


Matching < 100 


11.0 


14.1 


15.6 


13.4 


12.4 


10.4 


9.3 


Mt > 200 GeV 


8.6 


13.2 


15.2 


13.3 


12.4 


10.3 


9.3 


Charge survival efficiency (38%) 


3.3 


5.0 


5.8 


5.0 


4.7 


3.9 


3.5 



TABLE XII: Selection efficiencies (in %) for a single top squark in MC events. Initially, events with at least one isolated muon 
of good quality and pr > 20 GeV are required and a match between a muon and a central track is required based on x^- 



M(gaugino-like chargino) in GeV 


100 


150 


200 


250 


300 


Selection criteria 












Initial muon selection 


51.7 


54.7 


55.4 


54.3 


53.8 


Isolated muon 


41.3 


44.6 


45.7 


44.7 


44.9 


Trigger probability 


22.7 


24.3 


24.7 


24.2 


24.1 


Cosmic veto 


21.0 


20.8 


21.0 


19.9 


19.7 


LiL trigger gate, i^smt ^ o 


15.5 


13.8 


19 4 


10.0 


8.4 


PT > 60 GeV 


12.2 


13.1 


12.3 


10.0 


8.4 


(/5) < 1, Speed xVdof < 2 


10.4 


11.7 


11.0 


8.9 


7.4 


Matching x^ < 100 


10.3 


11.6 


10.9 


8.9 


7.4 


Mt > 200 GeV 


7.8 


10.6 


10.6 


8.8 


7.4 


TABLE XIII: Selection efficiencies (in %) 


for a single gaugino-like char; 


jino in simulated events 


. Initially, events with at least 


one isolated muon of good quality and pt 


> 20 GeV are 


required and a 


match between a muon 


and a central track 


is required 


based on x^- 












M(higgsino-like chargino) in GeV 


100 


150 


200 


250 


300 


Selection criteria 












Initial muon selection 


52.4 


55.6 


56.6 


56.1 


53.9 


Isolated muon 


42.2 


45.8 


47.0 


47.0 


45.3 


Trigger probability 


23.1 


25.0 


25.3 


25.3 


24.2 


Cosmic veto 


21.1 


23.4 


23.7 


21.3 


20.1 


LI trigger gate, Nsmt > 3 


15.9 


15.9 


14.5 


11.4 


9.1 


PT > 60 GeV 


12.7 


15.2 


14.4 


11.3 


9.1 


(/3) < 1, Speed xVdof < 2 


10.8 


13.4 


12.7 


10.1 


8.2 


Matching x^ <100 


10.8 


13.4 


12.7 


10.1 


8.1 


Mt > 200 GeV 


8.3 


12.2 


12.4 


10.0 


8.1 



TABLE XIV: Selection efficiencies (in %) for a single higgsino-like chargino in MC events. Initially, events with at least one 
isolated muon of good quality and pt > 20 GeV are required and a match between a muon and a central track is required 
based on x^- 
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masses 100-300 GeV. The BDT-output distributions for 
simulated stau lepton, top squark, and higgsino-like 
charginos, after being normalized to the expected num- 
ber of events, are shown in Figs. BTll^ in Appendix A. 
These BDT distributions are used as input to a modified 
frequentist limit calculation employing a log-likelihood 
ratio (LLR) test statistic [ill, \^ to obtain 95% C.L. 
cross section limits. To minimize the degrading effects 
of systematic uncertainties on the search sensitivity, the 
signal and background are fitted to the observed data by 
maximizing a likelihood function over the systematic un- 
certainties for both the background-only and the signal- 
plus-background hypotheses. 

C. Systematic Uncertainties 



FIG. 17: (color online) The transverse mass Mt for single 
muon events and higgsino-like chargino events with chargino 
masses 100 GeV, 200 GeV, and 300 GeV. The single muon 
event sample satisfies all of the selection criteria described in 
Table |X] except the Mt cut, which separates the data and 
background samples. All the selection criteria described in 
Table IXIVI have been applied to the charginos except the 
Mt > 200 GeV cut. The distributions have been normal- 
ized to the same number of events. 



Events containing muons with measured (/3) > 1 are used 
to define a signal-free region. Signal-free control events 
contain muons with {(3) > 1 and Mt < 200 GeV, and 
signal-free data contain muons with (/3) > 1 and Mt > 
200 GeV. If the number of background events is Nb, 
the number of signal- free control events is NsfCi and 
the number of signal- free data events is Nsfd, then the 
number of normalized background events, Nnbi can be 
expressed as 



NB 



Nf 



SFD 



NsFC 



(6) 



The key variables used for discrimination between sig- 
nal and background are (/?) and dE/dx. These vari- 
ables are anti-correlated for candidate tracks originat- 
ing from signal, but not for those originating from back- 
ground. Figure US] shows the adjusted d£'/da; as a func- 
tion of (/3) for simulated gaugino-like charginos, back- 
ground, and data. The variables (/?), speed significance, 
number of scintillation counter hits, dE/dx, dE/dx sig- 
nificance, and NsMT are used as inputs to a BDT. The 
BDT is trained to distinguish between signal and back- 
ground events using signal events from MC and back- 
ground events from data. Half the input events are used 
for training, while the other half are used as a test sam- 
ple to model the background and signal response of the 
BDT. The distributions of the BDT input variables are 
shown in Figs. [TH and HOI 

An example of the BDT-output distribution is shown 
in Fig. [5T] for simulated gaugino-like charginos for 



The sources of systematic uncertainties and their es- 
timation are similar to the search for CMLLP pairs, de- 
scribed in Sec. VI C. A parameter is varied within its un- 
certainty and the change is propagated through the entire 
analysis to produce a BDT distribution. The systematic 
uncertainties can be divided into two categories, "nor- 
malization" and "shape" uncertainties. Normalization 
uncertainties affect only the overall event rate whereas 
the shape uncertainties can also change the differential 
distribution. 

The systematic uncertainty due to uncertainty in lu- 
minosity (6.1%) and in muon reconstruction efficiency 
(2.1%) are normalization uncertainties. The other 
sources of normalization uncertainties are: the back- 
ground normalization uncertainty due to the choice of 
the cuts on (/3) (7.2%) and Mt (2.2%), the muon pT res- 
olution uncertainty (0.2%), the dE/dx correction uncer- 
tainty (< 0.1%), the dE/dx smearing uncertainty (0.2%), 
and the speed x^/dof correction uncertainty (0.4%). The 
systematic uncertainty for the speed x^/dof correction is 
determined by repeating the data to MC correction (de- 
scribed in Sec. lV A|) with a MC sample oiW fiv decays 
and then taking the difference between the two correc- 
tions as the uncertainty. The uncertainties due to the 
choice of PDF, and the dE/dx correction are the same 
as in Sec. VI C. The systematic uncertainties due to the 
width of the LI timing gate, and the timing simulation 
change with the output value of the BDT and therefore 
are shape systematics. The average uncertainty for the 
LI timing gate width is 4% and that for the timing simu- 
lation is 7%. These uncertainties are summarized in the 
third columns of Tables KVlIl and lXVIIII 



D. Results 

For the single CMLLP search with 5.2 fb~^ integrated 
luminosity, 95% C.L. upper limits on production cross 
sections for stau leptons, top squarks, gaugino-like, and 
higgsino-like charginos are shown in Table IXVI and in 
Fig. [211 Limits are set on production cross section of 
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FIG. 18: (color online) Adjusted dE/dx versus {j3) for (a) data events, (b) background (data events with Mt < 200 GeV), (c) 
MC sample for gaugino-like charginos with mass of 300 GeV. 



stau leptons from 0.05 pb to 0.006 pb, for stau lepton 
masses in the range between 100 GeV and 300 GeV. 
Pair-produced long-lived top squarks are excluded below 
masses of 285 GeV. If we only include the effects of initial 
hadronization of top squarks and do not include the ef- 
fects of charge flipping during their passage through the 
detector, the lower limit on long-lived top squark mass is 
found to be 305 GeV. Pair-produced long-lived gaugino- 
like charginos are excluded below masses of 267 GeV, and 
higgsino-like charginos below masses of 217 GeV. They 
are identical to those presented in 11|. Using the inter- 
section of the —1 (+1) standard deviation band with the 
NLO cross section, the mass limits shift down (up) by 1 
GeV for charginos and by 10 GeV for top squarks. 



VIII. RUN Ila RESULT 

The search for a pair of CMLLPs in 1.1 fb"^ of Run 
Ila integrated luminosity utilizes the TOF measurement 
in addition to other kinematic variables to select events 
with candidate CMLLPs [5i]. The CMLLP candidates 
in this analysis are staus, gaugino-like, and higgsino-like 
charginos. The principal background comes from mis- 
measured muons from the decays of the Z boson. Table 
IXVII presents the signal acceptance, the number of pre- 
dicted background events, and the number of observed 



events in this analysis. In the absence of any signal, limits 
of 206 GeV and 171 GeV are set on masses of gaugino-like 
charginos and higgsino-like charginos, respectively. Lim- 
its on the production cross section of stau leptons are set 
from 0.31 pb to 0.04 pb for the stau lepton mass range 
of 60 to 300 GeV. 



IX. COMBINATION OF RESULTS 



In the absence of observed signal in all three analyses, 
the searches for a pair of CMLLPs in 1.1 fb"^ of Run Ila, 
in 5.2 fb~^ of Run lib, and the search for a single CMLLP 
in 5.2 fb^^ of Run lib integrated luminosity, we combine 
the results to find limits on the production cross sections 
of stau leptons, gaugino-like charginos, and higgsino-like 
charginos. Due to the effect of hadronization and charge- 
flipping, the sensitivity of the search for single CMLLPs 
for top squarks is much better than the sensitivity of the 
search for a pair of CMLLPs. As a result, the top squark 
mass limits for the combination of the single CMLLP 
(with a veto for the events common with the pair search) 
and the pair CMLLP analyses do not show a significant 
improvement over the results from the single CMLLP 
analysis on its own and hence is not performed. 
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FIG. 19: (color online) Speed related input distributions to BDTs for the search for single CMLLPs. Background comes from 
events containing a muon with (/3) < 1 and Mt < 200 GeV. Signal is gaugino-like charginos of masses 100 GeV, 200 GeV, and 
300 GeV. The distributions are normalized to the same number of events. 



Method of Combination 



B. Systematic Uncertainties 



In the Run lib search for a pair of CMLLPs, the value 
of the requirement on the BDT output is optimized for 
each signal mass point. A similar procedure was used 
in the search for a pair of CMLLPs with Run Ila data. 
In the search for single CMLLPs in Run lib, the entire 
BDT distribution is used as input to a CLg limit set- 
ting method. To obtain results from the combination 
of the three analyses, the signal acceptance, the num- 
ber of predicted background events, and the number of 
observed events for each signal mass, for the two pair 
analyses (Tables IVlIll IXVH) along with the BDT distri- 
butions from the search for single CMLLPs (Figs.[5l] and 
[27ll29)) (see Appendix A) are used as inputs to the same 
CLs method. To avoid double counting of events, the 
datasets used for the two Run lib analyses are made sta- 
tistically independent by removing the events that have 
been selected for the search of CMLLP pairs from the 
dataset used to search for single CMLLPs. The different 
analyses are combined by summing the log-likelihood ra- 
tios over all the bins and all the analyses. 



Systematic uncertainties are treated as Gaussian dis- 
tributions and are applied to the expected number of 
signal and background events. Various sources of sys- 
tematic uncertainties along with their values for the four 
types of CMLLP signals that have been studied are listed 
in Table IXVIII and the systematic uncertainties on the 
background sample are listed in Table lXVlIIl All system- 
atic uncertainties except the luminosity uncertainty (ssj 
for the two searches for pairs of CMLLPs are treated as 
uncorrelated. The recent search for CMLLP pairs and 
the search for a single CMLLP are based on the same 
dataset. Therefore, the systematic uncertainties for these 
two analyses are correlated except for the uncertainties 
on background normalizations (the background samples 
are different in the two analyses) . The systematic uncer- 
tainties for the single CMLLP dataset after removal of 
the events containing CMLLP pairs are the same as those 
in the search for a single CMLLP analysis described in 
Sec. VII D. The shape systematic uncertainties however 
are updated after the removal of the common events. 
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FIG. 20: (color online) Distributions related to AE/Ax used in the BDT for the search for single CMLLPs. Background comes 
from events containing a muon with < 1 and Mt < 200 GeV. Signal is gaugino-like charginos of masses 100 GeV, 200 GeV, 
and 300 GeV. The distributions are normalized to the same number of events. 



C. Results 

Combined 95% C.L. cross section limits for stau lep- 
tons, gaugino-like, and higgsino-like charginos are shown 
in Table IXiXl and Fig. [23l Using the observed cross sec- 
tion and the theoretical NLO cross section, we set mass 
limits of 278 GeV for gaugino-like charginos and 244 GeV 
for higgsino-like charginos. Using the intersection of the 
— 1 (-1-1) standard deviation (cr) band on the NLO cross 
section shifts the mass limits down (up) by '^l GeV for 
the charginos. We do not have enough sensitivity to set 
a limit on the stau lepton mass and therefore we set an 
upper limit on production cross sections of stau leptons 
to be 0.04 pb to 0.008 pb for the stau lepton mass range 
of 100 to 300 GeV. 

X. SUMMARY 

A search for CMLLPs has been performed with the 
DO detector with 5.2 fb~^ integrated luminosity using 
two different strategies: a search for a pair of identified 
CMLLPs and a search for a single identified CMLLP in 
events expected to contain a pair of CMLLP's. These two 
searches are combined with the earlier search for CMLLP 
pairs with 1.1 fb~^ integrated luminosity. We use the 



central value of the theoretical cross section predictions 
to set 95% C.L. lower limits on the masses of top squarks 
and charginos and on the cross section of stau leptons. 

Using the combination of the three searches we set 
mass limits of 278 GeV for gaugino-like charginos and 
244 GeV for higgsino-like charginos. For stau leptons we 
set an upper limit of 0.04 pb to 0.008 pb on the pro- 
duction cross section for the mass range of 100 to 300 
GeV. 

In the search for single CMLLPs we exclude top 
squarks with masses below 285 GeV with a charge flip- 
ping probability of 38%. A combination of the analyses 
is not performed for the top squarks since improvement 
in the stop limit by combining the searches is negligible. 

Limits on the chargino cross sections are the most re- 
strictive limits to date, with about an order of magni- 
tude improvement over the previous DO result with 1.1 
fb^^ integrated luminosity The improvement in both 
the pair and the single CMLLP searches over the previ- 
ous results is due to the increased luminosity as well as 
the additional use of another key variable, the measured 
AE/Ax of the tracks. 
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FIG. 21: (color online) BDT-output distributions for simulated gaugino-like chargino masses 100-300 GeV in 50 GeV steps in 
the search for single CMLLPs. The distributions are normalized to the expected number of events. 



Acknowledgments 

We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
MON, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAFESF and FUNDUNESF (Brazil); DAE and DST (In- 



dia); Colciencias (Colombia); CONACyT (Mexico); NRF 
(Korea); FOM (The Netherlands); STFC and the Royal 
Society (United Kingdom); MSMT and GACR (Czech 
Republic); BMBF and DFG (Germany); SFI (Ireland); 
The Swedish Research Council (Sweden); and CAS and 
CNSF (China). 



[1] M. Byrne, C. Kolda, and P. Regan, Phys. Rev. D 66, 
075007 (2002). 

[2] K. Nakamura et al. (Particle Data Group), J. Phys. G 
37, 075021 (2010), see Big Bang Nucleosynthesis (rev.), 
Section 20.5 and references cited therein. 

[3] M. Eads, Ph.D. Thesis, Northern Illinois University, 
FERMILAB-THESIS-2005-24 (2004). 



[4] Y. Xie, Ph.D. Thesis, Brown University, FERMILAB- 

THESIS-2009-37 (2009). 
[5] V. Abazov et al. (DO Collaboration), Phys. Rev. Lett. 

102, 161802 (2009). 
[6] D. Acosta et al. (CDF Collaboration), Phys. Rev. Lett. 

90, 131801 (2003). 
[7] T. Aaltonen et al. (CDF Collaboration), Phys. Rev. Lett. 



28 



lOi 



={a) D0 5.2fb"' 



10-^ 



Observed Limit 
Expected Limit 
Expected ± 1 SD 
Expected ± 2 SD 
NLO Prediction 
NLO Uncertainty 




100 



150 200 250 
Stau Lepton Mass [GeV] 



300 



:{C) D0 5.2fb"' 



10-^ 



£1 

T 

■K- 10' 



10-' 



Observed Limit 
Expected Limit 
Expected ± 1 SD 
Expected ± 2 SD 
NLO Prediction 
NLO Uncertainty 




100 150 200 250 300 
Gaugino-Like Chargino Mass [GeV] 



-- (b) D0 5.2fb"' 



Observed Limit 
Expected Limit 
Expected + 1 SD 
Expected ± 2 SD 
Observed Limit, no CF 
Expected Limit, no CF 
NLO Prediction 
NLO Uncertainty 




100 150 200 250 300 350 
Top Squark Mass [GeV] 



400 



£1 

S 10-^ 



T 



10"' 



10-' 



D0 5.2 fb 



Observed Limit 
Expected Limit 
Expected ± 1 SD 
Expected ± 2 SD 
NLO Prediction 
NLO Uncertainty 




100 150 200 250 300 
Higgsino-Like Chargino Mass [GeV] 



FIG. 22: (color online) 95% C.L. limits on production cross sections of a pair of stau leptons, top squarks, gaugino-like charginos, 
and higgsino-like charginos, as a function of their masses from the search for one or more CMLLPs with Run lib data. ±1 SD 
and ±2 SD are the 1 and 2 standard deviation bands respectively around the expected limit curves. 
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TABLE XV: 95% C.L. cross section limits for stau lepton, top squark, gaugino-like, and higgsino-like charginos in the search 
for single CMLLPs in Run lib data. 
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Appendix A: BDT distributions 

The BDT-output distributions for stau, top squark, 
and higgsino-like charginos, after being normalized to the 
expected number of events, for the search of a pair of CM- 
LLPs are shown in Figs. [2414261 The BDT-output distri- 
butions for stau, top squark, and higgsino-like charginos, 
after being normalized to the expected number of events, 
for the search of a single CMLLP are shown in Figs. B71 - 

m 
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Mass (GeV) 


Signal Acceptance (%) 


Predicted Background 


Observed Data 


Stau lepton 








100 


5.56 ±0.11 ±0.41 


1.55 ± 0.49 ± 0.30 


1 


150 


12.3 ± 0.16 ± 1.27 


1.70 ±0.51 ±0.15 


1 


200 


13.9 ±0.17 ±1.11 


1.70 ±0.51 ±0.51 


1 


250 


13.3 ± 0.16 ± 1.25 


1.70 ±0.51 ±0.31 


1 


300 


11.7 ±0.15 ±1.34 


1.86 ±0.54 ±0.15 


2 


Gauginolike chargino 








100 


4.63 ±0.10 ±0.35 


1.55 ±0.49 ±0.31 


1 


150 


8.51 ±0.13 ±0.88 


1.24 ±0.44 ±0.11 


1 


200 


8.89 ±0.13 ±0.71 


1.86 ± 0.54 ± (< 0.01) 


1 


250 


7.40 ±0.12 ±0.70 


1.70 ±0.51 ±0.31 


1 


300 


5.88 ±0.11 ±0.68 


1.70 ±0.51 ±0.14 


2 


Higgsino-like chargino 








100 


4.94 ±0.10 ±0.37 


1.55 ± 0.49 ± 0.31 


1 


150 


8.91 ±0.13 ±0.92 


1.39 ±0.46 ±0.13 


1 


200 


9.56 ± 0.14 ± 0.76 


1.86 ± 0.54 ± (< 0.01) 


1 


250 


8.13 ±0.13 ±0.76 


1.70 ±0.51 ±0.31 


1 


300 


6.36 ±0.11 ±0.73 


1.70 ±0.51 ±0.14 


1 



TABLE XVI: The signal acceptance, number of predicted background events, and the number of observed events from the 
search for a pair of CMLLPs with 1.1 fb^^ of Run Ila integrated luminosity. The first error is statistical and the second one is 
systematic. 



Pair (1.1 fb~^) 



Pair (5.2 fb~^) 



Single (5.2 fb~^) 



Luminosity (±6.1%) 
Muon Reco. (±0.7%) 
PDF (±0.1-2.7%) 
Timing gate (±2.8-13%) 
Time simulation (±6-13%) 



Luminosity (±6.1%) 
Muon Reco. (±2.1%) 

PDF (±0.2%) 
Pt resolution (±2.8%) 
Timing gate (±2.4%) 
Time simulation (±2.8%) 

AE/Ax corr. (±0.1%) 
AE/Ax smearing (±0.6%) 
MDT Timing gate (±1.2%) 
Speed Asym. C'tjrr. ( = 1.06 lO.lVC) 



Luminosity (±6.1%) 
Muon Reco. (±2.1%) 

PDF (±0.2%) 
Pt resolution (±0.2%) 
Timing gate (shape) 
Time simulation (shape) 
AE/Ax corr. (±0.02%) 
AE/Ax smearing (±0.2%) 
Speed xVdof Corr. (±0.4%) 



TABLE XVII: Systematic uncertainties for signals for all three analyses. 



31 



Pair (1.1 fb"^) 



Pair (5.2 fb"^) 



Single (5.2 fb"^) 



Bkgd. Norm. (±9-28%) 



Luminosity (±6.1%) 
Muon Reco. (±2.1%) 
PDF (±0.3%) 
Pt resolution (±11.0%) 
Timing gate (±3.8%) 
Time simulation (±9.5%) 

dE/dx Corr. (±1.5%) 
dE/dx smearing (±4.9%) 
Bkgd. Norm. (±2.2%) 
Speed Asym. Corr. (±3.6%) 



dE/dx Corr. Uncertainty (±0.02) 
Bkgd. Norm.-(/3) (±7.2%) 
Bkgd. Norm.-MT (±2.2%) 



TABLE XVIII: Systematic uncertainties for background events for all three analyses. 



Mass (GeV) 



NLO cross section (pb) 



_o6s 
<^95 



(pb) 



(pb) 
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TABLE XIX: Combined 95% C.L. cross section limits for stau leptons, gaugino-like, and higgsino-like charginos from the three 
search strategies. 
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FIG. 23: (color online) Combined limits at 95% C.L. on production cross sections of a pair of stau leptons, gaugino-like 
charginos, and higgsino-like charginos as a function of their masses with Run Ila and Run lib data. ±1 SD and ±2 SD are the 
1 and 2 standard deviation bands respectively around the expected limit curves. 
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Output 



FIG. 24: (color online) BDT-output distributions for stau masses 100-300 GeV in 50 GeV steps for the search for a CMLLP 
pair with Run lib data. The distributions are normalized to the expected number of events. The selection requirement on the 
BDT value is shown with a green vertical line. 
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FIG. 25: (color online) BDT-output distributions for top squark masses 100-400 GeV in 50 GeV steps for the search for a 
CMLLP pair with the Run lib data. Distributions are normalized to the expected number of events. Selection requirement on 
the BDT value is shown with a green vertical line. 
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FIG. 26: (color online) BDT-output distributions for higgsino-like chargino masses 100-300 GeV in 50 GeV steps for the search 
for a CMLLP pair with the Run lib data. Distributions are normalized to the expected number of events. Selection requirement 
on the BDT is shown with a green vertical line. 



36 



10= r, ^ 
10^ r^^) 

10^ r 
^ 10^ r 

~~ 1 r 

B 10-' r 
10-^ - 

10-" r 

10-= r 
10-*^ |- 

-0.8 



D0 5.2 fb-^ 



• Data 

n Background 
-ffi- — TT, 100 GeV 



10' r, , 

10" 

10' r 
^ 10^ r 

a 10-' r 

10-" r 

10-= r 
10-*^ |- 

-0.8 



D0 5.2 fb-^ 



• Data 

n Background 
— TT, 200 GeV 




-0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 
BDT Output 



-0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 
BDT Output 



10= •- 

pb) D0 5.2fb-i 

io4 

,2 r 



1 r 
a 10-' r 
10-^ - 

10-" r 

10-= r 

-0.8 



• Data 

n Background 
l_ — Tp 150 GeV 



-0.6 -0.4 -0.2 -0 0.2 0.4 0.6 O.i 
BDT Output 



10= •- 

pd) D0 5.2fb-i 

104 
^ 10^ r 

1 r 



c 10-' r 
t 10-^ h 

10-" r 

10-= r 

-0.8 



• Data 

n Background 
— T7, 250 GeV 



-0.6 -0.4 -0.2 -0 0.2 0.4 0.6 O.S 
BDT Output 




• Data 

□ Background 
{ih — \, 300 GeV 



-0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 
BDT Output 



FIG. 27: (color online) BDT-output distributions for stau masses 100-300 GeV in 50 GeV steps in the search for single GMLLPs. 
The distributions are normalized to the expected number of events. 
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FIG. 28: (color online) BDT-output distributions for top squark masses 100-400 GeV in 50 GeV steps in the search for single 
CMLLPs. The distributions are normalized to the expected number of events. 
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FIG. 29: (color online) BDT-output distributions for higgsino-like chargino masses 100-300 GeV in 50 GeV steps in the search 
for single CMLLPs. The distributions are normalized to the expected number of events. 



